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Abstract 



^.^ ■ The Hamiltonian dynamics of classical planar Heisenberg model is numeri- 

(<— ^ ■ cally investigated in two and three dimensions. In three dimensions peculiar 

0^ . behaviors are found in the temperature dependence of the largest Lyapunov 

J>^. exponent and of other observables related to the geometrization of the dy- 

Q . namics. On the basis of a heuristic argument it is conjectured that the phase 

rS^ ■ transition might correspond to a change in the topology of the manifolds 



whose geodesies are the motions of the system. 
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On the basis of the ergodic hypothesis, Statistical Mechanics describes the physics of 
many-degrees of freedom systems by replacing time averages of the relevant observables 
with ensemble averages. In the present Letter, instead of using statistical ensembles, we 
investigate the Hamiltonian (microscopic) dynamics of a system undergoing a phase transi- 
tion. The reason for tackling dynamics is twofold. First, there are observables, like Lyapunov 
exponents, that are intrinsically dynamical. Second, the geometrization of Hamiltonian dy- 
namics in terms of Riemannian geometry provides new observables and, in general, a new 
interesting framework to investigate the phenomenon of phase transitions. 

The geometrical formulation of the dynamics of conservative systems |[l| was first used 
by Krylov in his studies on the dynamical foundations of statistical mechanics and sub- 
sequently became a standard tool to study abstract systems in ergodic theory. Several new 
contributions to this subject appeared in the last years 0-0. 

Let us briefly recall that the geometrization of the dynamics of A^-degrees-of-freedom 
systems deflned by a Lagrangian C = T — V, in which the kinetic energy is quadratic in the 
velocities: T = ^aijq'^q^ , stems from the fact that the natural motions are the extrema of 
the Hamiltonian action functional Sh = J Cdt, or of the Maupertuis' action Sm = 2 J T dt. 
In fact also the geodesies of Riemannian and pseudo- Riemannian manifolds are the extrema 
of a functional: the arc-length i = J ds, with ds"^ = Qijdq^dq^, hence a suitable choice of 
the metric tensor allows for the identiflcation of the arc-length with either Sh or Sm, and 
of the geodesies with the natural motions of the dynamical system. Starting from Sm the 
"mechanical manifold" is the accessible conflguration space endowed with the Jacobi metric 
{gj)ij = [E — V{{q})] ttij, where V{q) is the potential energy and E is the total energy. A 
description of the extrema of Hamilton's action Sh as geodesies of a "mechanical manifold" 
can be obtained using Eisenhart's metric ||^ on an enlarged conflguration spacetime {{q^ = 
t,q^, . . . , g^} plus one real coordinate g^^^), whose arc-length is 

ds^ = -2V(q) {dq^ + aijdq'dq^ + 2dq^dq^+^ . (1) 

The manifold has a Lorentzian structure and the dynamical trajectories are those geodesies 



satisfying the condition ds^ = Cdt"^, where C is a positive constant. In the geometrical 
framework, the (in)stabihty of the trajectories is the (in)stabihty of the geodesies, and it is 
completely determined by the curvature properties of the underlying manifold according to 
the Jacobi equation [Q 

^^+^,7c„.^^^-0, (2) 

whose solution J, usually called Jacobi or geodesic variation field, locally measures the 
distance between nearby geodesies; D/ds stands for the covariant derivative along a geodesic 
and R' jkm ^^^ ^^^ components of the Riemann curvature tensor. Using the Eisenhart metric 
(H) the relevant part of the Jacobi equation (H) is 0,^ 

^ + i?Vo'/' = 0, z = l,...,iV (3) 

where the only non- vanishing components of the curvature tensor are -RojOj = d'^V/dqidqj. 
Equation (^ is the tangent dynamics equation which is commonly used to measure Lyapunov 
exponents in standard Hamiltonian systems. Having recognized its geometric origin, in ref. 
1^ we have devised a geometric reasoning to derive from Eq. (^) an effective scalar stability 
equation that independently of the knowledge of dynamical trajectories provides an average 
measure of their degree of instability. This is based on two main assumptions: i) that 
the ambient manifold is almost isotropic, i.e. the components of the curvature tensor - 
that for an isotropic manifold (i.e. of constant curvature) are Rijkm = i^oigikOjm ~ gimQjk)-, 
kq = const - can be approximated by Rijkm ~ i^if){.gik9jm — QimQik) along a generic geodesic 
7(t); ii) that in the large A^ limit the "effective curvature" n{t) can be modeled by a gaussian 
and ^-correlated stochastic process. The mean kq and variance o"k of n{t) are given by the 
average and the r.m.s. fluctuation of the Ricci curvature kR = Kr/N along a geodesic: 
kq = {Kr)/N, and a^ = {{Kr — {Kr)Y) /N respectively. The Ricci curvature along a 
geodesic is defined as Kr = Rij^-§f/i-jf-^), where Rij = R'^ikj is the Ricci tensor; in the 
case of Eisenhart metric it is Kr = AV = Yli=i d'^V/dqf. The final result is the replacement 
of Eq.(^ with the aforementioned effective stability equation which is independent of the 
dynamics and is in the form of a stochastic oscillator equation P,|^ 



^ + ^W^ = 0, (4) 

where ip"^ oc |Jp. The exponential growth rate A of the solutions of Eq. @, which is 
therefore an estimate of the largest Lyapunov exponent, can be computed exactly: 



^4-^' ^H-^-Vl?^--^- I (^) 



where r = Ti^fKQ/{2^ kq{kq + a^) + ikJi^)] in the limit (Ti^/ kq <^ 1 one finds A oc cr^. Details 
can be found in Refs. f^J^. 

In our geometric picture chaos is mainly originated by the parametric instability 
activated by the fluctuating curvature "felt" by the geodesies. On the other hand, the 
average curvature properties are statistical quantities like thermodynamic observables. This 
means that there exists a non-trivial relationship between dynamical properties (Lyapunov 
exponents) and suitable static observables. Generic thermodynamic observables have a 
non-analytic behaviour as the system undergoes a phase transition. Hence the following 
question arises naturally: "Is there any peculiarity in the geometric properties associated 
with the dynamics, and thus in the chaotic dynamics itself, of systems which exhibit an 
equilibrium phase transition?" And in particular, do the curvature fluctuations and/or 
the Lyapunov exponent show any remarkable behaviour in correspondence with the phase 
transition itself? We address this question considering a system of planar classical "spins" 
(rotators) Sj = {co^ (pi,sunpi) defined on a d- dimensional lattice U^. The Hamiltonian is 

n{{^.^}) = \Y,^l + v{{^,}), (6) 

i 

where (pi and ttj are the canonically conjugated angle and angular momentum of the "spin" 
on the i-th lattice site. The interaction is given by {{ij) stands for nearest-neighbour sites) 

V = - Y. (cos(^,-^,)-l) , (7) 

which is the Heisenberg XY potential. We consider c? = 2, 3. The potential (|^ is invariant 
under the action of the continuous group 0(2), hence — in the limit N —* oo — we expect 



a second order phase transition only in d = 3 and a Kosterlitz-Thouless transition in d = 2. 
The equations of motion derived from the Hamiltonian (O) have been numerically integrated 



using a symplectic algorithm |T^, with random initial conditions at equipartition (energy 
equally shared among the degrees of freedom) and at several values of the energy density 
e = E/N. At each e we measured the corresponding temperature T as the time average 
of the kinetic energy per degree of freedom. The temperature behavior of internal energy, 
specific heat and vorticity, computed as time averages instead of ensemble averages, led us 
to estimate a critical temperature Tc ^ 0.95 in d = 2, in agreement with the already existing 
estimates ^^, and Tc ~ 2.15 in the d = 3 case. In Figs. |I] and |^ the values of the largest 



Lyapunov exponent, numerically computed using the standard algorithm |T2[, are plotted vs. 
the temperature T and are compared to their corresponding analytic estimates obtained by 
means of Eq. (^, where kq and a^ are computed as time averages. The agreement between 
theoretical predictions and numerical data is very good; in an intermediate temperature 
range a "renormalization" of Kq is necessary, as already discussed in Ref. for the one- 
dimensional case. 

In the d = 2 case Ai(T) displays a rather smooth pattern in the transition region (see 
the inset of Fig. 1), whereas in the d = 3 case, at T ~ 2.15 = T^., the behavior of Ai(T) 
clearly shows a neat departure from its intermediate regime of linear growth, as can be seen 
in the inset of Fig. ^ where the transition region is magnified and linear scales are used. No 
evidence of a possible divergence of Ai(T) is found as T ^ T^, at variance with the results 
reported in Ref. |]13[, though a very different model is considered therein. In this respect our 



results for Ai(T) are closer to those found for a liquid-solid first-order transition |T^ and for 
other models 111511 . 



Let us now turn to the hidden geometry of the dynamics and in particular to the complex 
landscape of the ambient manifold whose deviation from isotropy — quantified by a^ — is 
directly responsible for dynamical chaos. The comparison of Figs, ^and^, where ko{T) and 
(Jk{T) are reported ioT d = 2 and d = 3 respectively, evidences a remarkable feature of the 
curvature fluctuations: a singular (cusp-like) behavior of (Jk{T) shows up in correspondence 



with the second order phase transition and <Jk(T) is sharply peaked at Tc, whereas in absence 
of symmetry breaking {d = 2) no singular behaviour of ct^{T) is present. This behavior of the 
curvature fluctuations is very intriguing. In fact a singular behavior of the curvature fluctu- 
ations can be reproduced in abstract geometric models which undergo a transition between 
different topologies at a critical value of a parameter that can be varied continuously. Let us 
consider for instance the families of surfaces of revolution immersed in M? defined as follows: 
^e = (/elw) cos f ,/£(«) sin !>,«), where m, f are local coordinates on the surface {v E [0,2tt] 
and u belongs to the domain of definition of fs), fe{u) = iV^ + v? — u"^ , e G [^min? +00), 
and emin = —\- There is a critical value of the parameter, £ = ^c = 0, corresponding to 
a change in the topology of the surfaces. In particular the manifolds JF^ are diffeomorphic 
to a torus T^ when e < and to a sphere §^ when e > 0. Computing the Euler-Poincare 
characteristic x by means of the Gauss-Bonnet theorem |jT6|, one finds xi^s) = if e < 0, 
and x(-^e) = 2 otherwise. 

Let M be a generic member of the family JF^, and let us define the fluctuations of the 
gaussian curvature K (see e.g. Ref. [p6| for the definition of K) as a"^ = (K^) — {K)"^ = 



A ^ J K"^ dS — A ^{J K dS)"^ where A is the area of M and dS is the invariant surface 

M M 

element. This family of surfaces exhibits a singular behaviour in the curvature fluctuation 
0" as e — i> 6c, as shown in Fig. ^. This is remarkably similar to the cusp-like behavior of the 



Ricci curvature fluctuations cr«;(T) of the XY model in rf = 3 that are peaked at Tc |17 . 
At heuristic level, these results suggest that a phase transition might correspond to a major 
topology change in the manifolds underlying the motion. We conjecture that the family of 
"mechanical manifolds" (each one being in one-to-one correspondence with a value of T) 
splits, at Tc, into two subfamilies of manifolds that are not diffeomorphic (being perhaps of 
different cohomology type). 

The relevance of topological concepts for the theory of phase transitions has been al- 



ready rigorously demonstrated in a rather abstract context (see Ref. [0); the present work 
suggests that also topological properties of the manifolds underlying the microscopic (Hamil- 
tonian) dynamics could be relevant to second order phase transitions. 
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FIGURES 
FIG. 1. Lyapunov exponent Ai vs. energy density e for the d = 2 case. Numerical results 

correspond to lattice size: A*" = 10^ (starred open squares), N = 20^ (open triangles), N = 40^ 

(open stars), N = 50^ (open squares), N = 100^ (open circles). Full squares are analytic results 

according to Eq.(4); dots are analytic results without correction (see text). In the inset symbols 

have the same meaning. 

FIG. 2. Lyapunov exponent Ai vs. e for the d = 2> case. Numerical results with lattice size: 
N = 10^ (open squares), N = 15^ (open star). Analytic results are represented by full circles; dots 
are analytic results without the correction mentioned in the text. Inset: N = 10^ (open squares), 
N = 15^ (open circles). 

FIG. 3. d = 2 case. Time averages, at iV = 40^, of Ricci curvature (open circles) and its 
fluctuations (full circles). Solid lines are analytic estimates obtained from a high temperature 
expansion. 

FIG. 4. d = 3 case. Time averages, at iV = 10^, of Ricci curvature (open triangles) and its 
fluctuations (full triangles). Open circles and full rhombs refer to a lattice size of A^ = 15^. Solid 
lines are microcanonical analytic estimates obtained from a high temperature expansion. The 
appearence of a cusp-like behavior of curvature fluctuations is well evident at £c- 

FIG. 5. Fluctuations amplitude, o", of Gauss curvature of a family of surfaces parametrized by 
e. For graphical reasons e is shifted by its minimum value \emin\ = 0.25, thus the cusp corresponds 
to e = 0, the critical value separating two families of different Euler characteristic x i-6- of different 
topology. 
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